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Abstract. 

We present numerical simulations of circular Couette flow in axisymmetric and 
fully three-dimensional geometry of a cylindrical annulus inspired by Princeton MRI 
liquid gallium experiment. The incompressible Navier-Stokes equations are solved 
with the spectral element code Nek5000 incorporating realistic horizontal boundary 
conditions of differentially rotating rings. We investigate the effect of changing rotation 
rates (Reynolds number) and of the horizontal boundary conditions on flow structure, 
Ekman circulation and associated transport of angular momentum through the onset 
of unsteadiness and three-dimensionality. A mechanism for the explanation of the 
dependence of the Ekman flows and circulation on horizontal boundary conditions is 
proposed. 
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1. Introduction 

The phenomenon of Ekman circulation (EC) occurs in most if not all rotating flows 
with stressed boundaries that are not parallel to the axis of rotation. The manifestation 
of EC ranges from wind-driven ocean currents (Batchelor 1967), to the accumulation 
of the tea leaves at the bottom of a stirred cup (see, e.g., Alpher & Herman 1960). 
One of consequences of EC and of the associated Ekman flows is greatly to to enhance 
mixing and transport and in particular, the transport of angular momentum, above 
the values due to viscosity alone. Traditionally, Ekman flows are explained in terms 
of action of Coriolis forces in the Ekman layers along the rotating stressed boundaries 
(Greenspan 1968). 

There are circumstances when the presence of EC has undesirable effects. For 
example, this is the case in laboratory experiments to study the development of 
magneto-rotational instability (MRI) in liquid metals (see a monograph edited by Rosner 
et al. 2004). The MRI instability is important in astrophysics where it is believed to 
lead to turbulence in magnetized accretion disks (Balbus 2003). Many of the features 
of the MRI and its associated enhancement of angular momentum transport (AMT) 
can be studied experimentally in magnetized flows between rotating coaxial cylinders. 
In these experiments, the rotation rates of the cylinders are chosen in such a way that 
the fluid's angular momentum increases outwards so that the resulting rotational profile 
is stable to axisymmetric perturbations (so-called centrifugally stable regime). The 
presence of a weak magnetic field can destabilize the basic flow, provided the angular 
velocity increases inward, and lead to an enhancement of outward AMT. 

In an ideal situation, the basic state consists of circular Couette flow (CCF), and 
the outward transport of angular momentum in the absence of magnetic fields is solely 
due to viscous effects. The presence of a magnetic field would destabilize the basic 
flow through the effects of MRI and lead to a measurable increase of AMT. In practice, 
this ideal case can never be realized in laboratory experiments because of horizontal 
boundaries. The presence of these boundaries drives an EC that enhances AMT even 
in the absence of magnetic effects. In order to study the enhancement of AMT due to 
MRI it is crucial to be able to distinguish the effects that are magnetic in origin from 
those that are due to the EC. One possibility is to make the cylinders very tall so the 
horizontal boundaries are far removed from the central region. This approach, however, 
is not practical owing to the high price of liquid metals. 

The alternative approach is to device boundaries in such a way that the resulting 
EC can be controlled and possibly reduced. For example, attaching the horizontal 
boundaries to the inner or outer cylinder results in dramatically different flow patterns. 
Another possibility could be to have the horizontal boundaries rotating independently 
of inner and outer cylinder. Goodman, Ji and coworkers (Kageyama et al. 2004, Burin 
et al. 2006, Ji et al. 2006) have proposed to split the horizontal boundaries into two 
independently rotating rings whose rotational speeds are chosen so as to minimize the 
disruption to the basic CCF by secondary Ekman circulations. Indeed this approach has 
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been implemented in the Princeton's MRI liquid gallium experiment (Schartman 2008). 
In any case, no matter how the horizontal boundaries are implemented it is important to 
understand what kind of EC patterns arise before the magnetic effects are introduced. 

In the present paper we address this issue by studying the effects of horizontal 
boundary conditions on CCF numerically. We study both axisymmetric and fully 
three-dimensional geometries and investigate the effects of changing rotation rates 
(Reynolds number) through the onset of unsteadiness and three- dimensionality. The 
next section (section 2) describes the formulation of the problem and gives an account 
of numerical aspects of its solution technique including a brief description of the 
spectral element code Nek5000 (Fischer et al. 2008). The section 3 starts with an 
explanation of flow behaviour due to horizontal boundary conditions, i.e. CCF, Ekman 
and disrupted Ekman circulation due to periodic horizontal boundaries, 'lids' and 'rings', 
correspondingly (section 3.1). Then the paper proceeds with description of comparison 
of our results with the experimental data (section 3.2) followed by an examination of 
torque and AMT (section 3.3). Finally, we draw conclusions and describe future work 
in section 4. 



2. Problem Formulation and Numerical Method 



2.1. Formulation 



We study the flow of an incompressible fluid with finite (constant) kinematic viscosity 
v in a cylindrical annulus bounded by coaxial cylinders. The cylinders have the radii 
R\ and R* 2 (-R* < R 2 ) and rotate with angular velocities Q* and Q 2 , respectively. The 
annulus is confined in the vertical direction by horizontal boundaries at distance H* 
apart. The formulation of the problem in cylindrical coordinates (r, 6, z) with the scales 
for characteristic length L and velocity U, 



L 



K 2 



R{ U = Q{Rl - Q*R* (1) 

and therefore, with the relationship between dimensional variables (with asterisk) and 
non-dimensional radius, height, velocity vector V, time and pressure given by 



[r*,z*,V*,t*,p* 



Lr,Lz,UV,-t,pU 2 p 
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correspondingly, results in the following non-dimensional incompressible Navier-Stokes 
equations: 
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where p is a constant fluid density and Reynolds number Re is defined as 

Re = UL = (niRi-n^jR^-Rt) 

V V 

while the scalar advection operator due to a vector field V and laplacian of a scalar 
function S(r, z) are given by 

v)s _ v ds_ i^as y ds_ AS _<Ps_ ldS_ l_c?S_ c?S_ 

r dr r 89 z dz dr 2 r dr r 2 d6 2 dz 2 

The initial conditions for the flow in the annulus and boundary conditions at the 
cylinder surfaces r = Ri and r = R 2 are 

V r = V z = V e = rn(r) (9) 

where non-dimensional angular velocity f2(r) is given by circular Couette flow (CCF) 
profile 

^- A + 7 2 A ~ Rl-R 2 B ~ R\-R\ (10) 
At the horizontal boundaries z = and z = H, two types of the boundary conditions 
have been considered, namely, lids and rings, given by (9) where angular velocity Q(r) 
is equal to 

r = Ri 
i?i < r < R l2 




R 12 <r<R 2 ^ 



r = R 2 

Here i? i2 is the radial location of the boundary between the inner and outer rings, and 
and are angular velocities of inner and outer rings, correspondingly. Inspired 
by Princeton MRI liquid gallium experiment (Schartman 2008), the non-dimensional 
angular velocities and cylinder height as well as cylinder and ring boundary radii used 
in this study are given in table 1 in addition to the dimensional parameters involved 
in comparison with the experiment (subsection 3.2). In the cases with lids, angular 
velocities SI3 and are equal to the angular velocity of the outer cylinder Q 2 while in 
the cases with rings they turn out to be close to the values of CCF profile (10) taken at 
the middle of radii of the corresponding rings. 



2.2. Numerical Technique 

The axisymmetric version of equations (3-6) and fully three-dimensional version in 
cartesian coordinates has been solved numerically with the spectral-element code 
Nek5000 developed and supported by Paul Fischer and collaborators (see Fischer 
et al. 2008, Fischer et al. 2007, and references within). 

The temporal discretization in Nek5000 is based on a semi-implicit formulation 
in which the nonlinear terms are treated explicitly in time and all remaining linear 
terms are treated implicitly. In particular, we used either a combination of /cth- 
order backward difference formula (BDFA;) for the diffusive/so lenoidal terms with 
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26 




Table 1: The geometry and rotation parameters for the computational cases with lids 
and rings at Re = 6190 and experimental setup with lids at Re = 9270 along with the 
drawing of the cut of 3D computational mesh at 9 = for the case with rings. Note 
the clustering of the gridlines at the boundaries of the spectral elements whose location 
and dimensions are chosen to resolve efficiently boundary layers and 'step' changes in 
angular velocity between cylinders and rings. 



extrapolation (EXTA; — 1) for the nonlinear terms or the operator-integration factor 
scheme (OIFS) method where BDF/c is applied to the material derivative with the 
explicit fourth-order Runge-Kutta scheme being used for the resulting pure advection 
initial value problem. 

With either the BDF/c/EXT/c — 1 or OIFS approach, the remaining linear portion 
of time advancement amounts to solving an unsteady Stokes problem. This problem 
is first discretized spatially using spectral-element method (SEM) and then split into 
independent subproblems for the velocity and pressure in weak variational form. The 
computational domain is decomposed into K non-overlapping subdomains or elements, 
and within each element, unknown velocity and pressure are represented as the tensor- 
product cardinal Lagrange polynomials of the order N and N — 2, correspondingly, 
based at the Gauss-Lobatto-Legendre (GLL) and Gauss-Legendre (GL) points. This 
velocity-pressure splitting and GLL-GL grid discretization requires boundary condition 
only for velocity field and avoids an ambiguity with the pressure boundary conditions 
in accordance with continuous problem statement. 

The discretized Stokes problem for the velocity update gives a linear system which 
is a discrete Helmholtz operator. It comprises the diagonal spectral element mass matrix 
with spectral element Laplacian being strongly diagonally dominant for small timesteps, 
and therefore, Jacobi (diagonally) preconditioned conjugate gradient iteration is readily 
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employed. Then the projection of the resulting trial viscous update on divergence-free 
solution space enforces the incompressibility constraint as the discrete pressure Poisson 
equation is solved by conjugate gradient iteration preconditioned by either the two- level 
additive Schwarz method or hybrid Schwarz/multigrid methods. Note that we used 
dealising/overintegration where the oversampling of polynomial order by a factor of 
3/2 was made for the exact evaluation of quadrature of inner products for non-linear 
(advective) terms. 

The typical axisymmetric case with rings at high Reynolds number of Re = 6200 
(see figure 4b) required the spacial resolution with polynomial order N = 10 and number 
of spectral elements K = 320 (cf. drawing for table 1) and was computed with timestep 
At = 10~ 3 for the duration of t ~ 300, while the axisymmetric run with lids at the same 
Re (figure 4a) had N = 8, K = 476, At = 5 x 10~ 3 and t ~ 500. The corresponding 
three-dimensional cases with rings and lids had N — 11, K — 9600, At = 6.25 x 10 -4 , 
t ~ 280 and N = 9, K = 14280, At = 6.25 x 10~ 4 , t ~ 180, respectively. Note that in 
order to facilitate time advancement and minimize CPU requirements, the final output 
from another cases, e.g. with lower Reynolds number Re, was used as initial conditions 
for some of the computations with higher Re, and the corresponding axisymmetric 
cases with small random non-axisymmetric perturbation was a starting point for most 
of our fully 3D computations. Apart from CPU savings, the usage of the perturbed 
axisymmetric solution obtained in cylindrical formulation (3-6) as initial condition for 
3D computations at low Reynolds numbers (Re = 620) served as an additional validation 
of the code setup due to the convergence of the fully 3D results computed in cartesian 
formulation back to the unperturbed axisymmetric steady state initial condition (see 
also subsection 3.3). 

Finally, the step change of angular velocities that mimics its transition in the gaps 
or grooves between the cylinders and horizontal boundaries as well as between the 
inner and outer ring in Princeton MRI liquid gallium experiment (Schartman 2008) was 
modelled within one spectral element of the radial size L g = 0.020 by ramping power 
law function of radius with an exponent that was varied in the range from 4 to iV — 1 
without noticeable effect on the flow. 

3. Results 

Let us first start with examination of effects of horizontal boundary conditions on flow 
pattern in general and Ekman circulation in particular before moving to a comparison 
with the experiment and examination of angular momentum transport in the cylindrical 
annulus. 

3.1. Horizontal Boundary Effects 

Here we contrast two type of horizontal boundary conditions with an ideal baseline case 
of circular Couette flow (CCF). Being zero in the ideal case, we argue that the unbalance 
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Figure 1: The CCF azimuthal velocity Vg (- -), angular velocity Qc {~ ~~ ) ; axial 

angular momentum C = rVg ( — ■ — ) and minus pressure Pq ( ) versus radius. 

Note the monotonically increasing angular momentum and decreasing angular velocity 
with radius for centrifugally stable circular Couette flow where the 'centrifugal' rotation 
balances the centripetal pressure gradient leading to zero radial and axial velocities. 



between 'centrifugal' rotation and centripetal pressure gradient determines the fate of 
the radial flow along horizontal boundaries in the cylindrical annulus. 

In the ideal case of CCF, the sheared circular motion is balanced by centripetal 
pressure gradient. To be precise, the ideal CCF is the following exact solution of 
equations (3-6) for periodic (or stress-free) horizontal boundary conditions: 

V r = V z = V e = rn c (r) = Ar + - 

r 

rr y2 J^ r 1 g2 

Voir) = / -^dr = -— - 7r = + 2AB log r + Const (12) 
J r 2 2r z 

Here the constant A given by equation (10) is proportional to the increase in axial 
angular momentum, 

C = rV e = Q r 2 (13) 

outward between the cylinders while the constant B is set by shear-generating angular 
velocity drop between them. The figure 1 shows CCF azimuthal velocity Vg (dashed), 
angular velocity Qc (solid), axial angular momentum C (dash-dotted) and negative of 
pressure, —pc (dotted,) for the non-dimensional parameters given in table 1. Since we 
are primarily interested in further MRI studies, the baseline flow has to be centrifugally 
stable, i.e. with angular momentum C increasing outward (for Q > 0), and therefore, 
satisfying Rayleigh criterion 

dC 2 



dr 



> (14) 
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which is the case in this ideal CCF (dash-dotted line in figure 1). In order to maintain 
rotation with shear in this virtual experiment with periodic horizontal boundaries, the 
positive axial torque Tc 



T c = f (fx(dA-T)) = [ [ dzdO^ 

Ja^ ' z Jo J0 vi i Ve=rUc 

has to be applied to the inner cylinder while the outer cylinder is kept from shear-free 
solid body rotation (£l(r) = Q± = A, B = T = 0) by negative torque — 7c- Note that in 
the above equation (15), r is the non-dimensional shear stress tensor (see also Appendix 
B). 

3.1.1. Ekman Circulation with 'Lids' In practice, the ideal CCF can never be realized 
in laboratory experiments because of horizontal boundaries. The simplest realizable 
configuration is the one we refer to as 'lids' when horizontal boundaries are coupled to 
the outer cylinder (Q 3 = fl 4 = f2 2 ). To see how flow changes in the presence of lids that 
rotate with outer cylinder, let us imagine that these lids were inserted impulsively into 
fluid with ideal CCF profile given by equation (12) and plotted as a solid line in figure 1 
for Qi and Q 2 from table 1. Keeping the most important terms in the axisymmetric 
form of the equation (3) gives 

dV r 2 dp 1 d 2 V r 

-W = nr -^ + Te^ + --- (16) 
where we used Vq = rfl. For the initial condition of CCF (12), the left-hand side of 
equation (16) is equal to zero everywhere outside the lids which is also consistent with 
zero radial flow V r = 0. This zero radial flow also results in zero diffusion term j^^r in 
equation (16) and zero net radial force Q 2 r — The latter results from the exact CCF 
balance between (positive) 'centrifugal' rotation term fl c r and (negative) centripetal 
pressure gradient term — |^ in equation (16). 

Instead of initial ideal CCF angular velocity Qc (12), the flow next to the lids now 
rotates with a smaller angular velocity of the outer cylinder (Q 2 — ^3 — ^4 < ^c)- 
However, the centripetal pressure gradient is still set by the bulk rotation of the rest of 
the fluid and therefore, becomes suddenly larger than the 'centrifugal' rotation of fluid 
next to the lids, i.e. |^ = VL c r > Vl 2 r. As a result of this angular momentum deficit 
of near- wall fluid, the centripetal pressure gradient prevails over rotation term in (16). 
Therefore, the net radial force becomes non-zero and negative, Q 2 r — |^ < 0, resulting 
in negative sign of ^ (16) and therefore, in a formation of the Ekman layer with an 
inward radial flow (V r < 0) in the vicinity of the lids. 

The figure 2(a) confirms that the net radial force near, e.g. the lower horizontal 
surface z — 0, Vt\r — §f | 2 _ (dashed) is negative, as well as the scaled poloidal wall 
shear -== ^¥ L \ n . The latter means that the z-derivative of V r is negative at the lower 
lid which in turn results in a decrease of radial velocity with the increase of height z 
from noslip zero value at the lid, V r = (9) to negative values associated with the 

z=0 

inward Ekman flow. Thus the deficit of angular momentum in the near-wall fluid of the 
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(b) Rings 

Figure 2: Steady state scaled radial wall shear (— — ) and near- wall net radial force 
(- -) for Re = 620 in the case of lids (a) and rings (b). The definite negative net 
radial force in the lids case (a) result in the inward radial Ekman flow with negative 
radial wall shear being disrupted in the case of rings (b) by the alternating sign of the 
net radial force that correlates well with the sign of the wall shear and thus with the 
alternating directions of Ekman flows. 
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Ekman layer results in unbalanced centripetal pressure gradient set by the bulk rotation 
of the rest of the flow outside the layer and drives the Ekman flow radially inward. 

To summarize, the presence of slower rotating lids disrupt the initial ideal CCF 
equilibrium between centrifugal rotation and centripetal pressure gradient set by 
rotation of, respectively, lids and bulk of the flow. This leads to the negative net radial 
force, f2 2 r — < and inward Ekman flow, V r < owning to ^ < 0. As time grows, so 
does the magnitude of negative radial velocity in the Ekman layer and, eventually, the 
diffusion term ^^pr (16) in the Ekman boundary layer of the width Az ~ 0(VR~e) 
becomes of the same order (i.e. ~0(1)) as the net radial force that results from two 
other terms in (16). Thus the diffusion effects finally balance the rotation momentum 
deficit of the fluid in Ekman boundary layers near the lids in the saturation steady state 
(see also Appendix A). 

To check consistency of this argument, the saturation magnitude of ^ across the 
layer is verified to be more than an order of magnitude smaller than the corresponding 
I 2 which confirms that saturation centripetal pressure gradient —I 2 is indeed set by 
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Figure 3: Azimuthal velocity Vq versus radius for the circular Couette flow (— — ) and 
instantaneous axisymmetric profiles at z = -j in the case of lids for the series of Reynolds 

numbers Re = 620 ( ), 1900 (dash-triple-dot) and 6200 ( — • • — ), and in the case 

of rings for the same Reynolds numbers: ( ), ( — • — ) and (- -), respectively. 

The Ekman- circulation induced momentum deficiency in azimuthal velocity profiles in 
the cases with lids is greatly diminished by the particular choice of angular velocities of 
independently rotating rings. 



the bulk rotation of the fluid outside the Ekman boundary layers. The saturation 
bulk rotation can be illustrated by the instantaneous saturation profiles of azimuthal 
velocity shown in figure 3 for the steady cases with lids for Re = 620 (dashed) and 
1800 (dash-triple-dot), and unsteady case with lids of Re = 6200 (dash-double-dot) 
at z = j. Interesting that for the range of Reynolds numbers considered, the effect 
of the increase of Reynolds number is minor in comparison with significant azimuthal 
momentum deficiency resulted from the change of horizontal boundary conditions from 
initial ideal CCF (solid) to the cases of Ekman flows over lids. 

Owing to momentum deficiency of the near wall fluid, the higher centripetal 
pressure gradient drives the inward Ekman flows along the lids that result in EC in 
the cylindrical annulus. In order to further illustrate the phenomenon of EC due to 
horizontal boundaries, we have plotted the contours of azimuthal vorticity ujq and vector 
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(a) Lids 




(b) Rings 



Figure 4: Instantaneous contours of azimuthal vorticity and vector field of poloidal 
velocity for Re = 6200 in the case of lids (a) and rings (b). The Ekman circulation and 
outward radial jet near the midplane in the case of lids is severly disrupted in the setup 
with rings due to alternating inward- outward Ekman flows. 



plot of poloidal velocity (V r , V z ) in figure 4(a) for the case of Re 
given by 



6200 with the former 
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dV z 
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where noslip conditions (9) along the walls has been used. Here the vorticity contours 
are coloured from blue (u$ < 0) to red (cjg > 0). Note that this change of colours from 
blue to red (through green colour whenever it is visible) shows the locus of zero vorticity 
that gives approximate location of a jet or jet-like features in the flows along the lids at 
z = or z = H with extremum in V r and in the flows along the cylinders at r — Rx or 
r = i?2 with minimum or maximum in V z (17). In figure 4(a), we observe the Ekman 
boundary layers along the lids with vorticity contours changing their colours from blue 
{ojq < 0) to red at the lower lid and from red (ujq > 0) to blue at the upper lid. In both 
instances, this change of colour shows the locus of minimum in V r < or location where 
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inward Ekman flow is the strongest. Similarly, the change of colours near the inner 
cylinder surface shows the opposing vertical jet-like flows along the inner cylinder that 
merge near the midplane z = H/2 and owing to continuity (6), form strong outward 
radial jet. At these high Reynolds numbers, beyond Re ~ 1800, the radial jet becomes 
unsteady and starts to oscillates breaking into pairs of vortices or to be precise, into 
pairs of vortex rings that move toward the lids drawn by the mass loss in the Ekman 
layers and thus closing the EC cycle. 

Summing up the flow pattern in the case of lids, we conclude that because of the 
deficit of rotation momentum in Ekman layers, the fluid is pushed centripetally in these 
layers along the lids and further along the inner cylinder with the subsequent formation 
of the strong outward radial jet that eventually transport fluid back to the lids and 
closes the cycle of the EC (see also Appendix A and section 3.3). 

3.1.2. Ekman Circulation Disruption Due to 'Rings' When each horizontal boundary 
is split into a pair of rings that rotate independently with the angular velocities Q 3 and 
^4 (table 1), the bulk rotation and resulting centripetal pressure gradient is restored 
back to that of the CCF. The restored profiles of azimuthal velocity in the cases with 
rings are shown in figure 3 for the same Reynolds numbers as for the cases with lids, 
namely for the steady cases of Re = 620 (dotted) and 1800 (dash-dot) and unsteady 
case of Re = 6200 (long dash). Along with the restoration of the bulk rotation back to 
that of the CCF, we observe other major differences between the cases with lids (a) and 
rings (b) in the flow field structure illustrated in figure 4. 

Instead of a single outward radial jet and inward Ekman flows along the lids, 
figure 4(b) shows alternating inward-outward Ekman flows along the rings that, as we 
describe below, produce strong vertical jets near r = R\ 2 and a weaker outward radial 
jet near midplane z = H/2. The alternating inward-outward Ekman flows along, e.g., 
the lower inner and outer rings (z — 0), are also evident in figure 2(b) where the scaled 
poloidal shear ^== ^f| z _ (solid) is plotted as a function of radius for Re = 620. The 
radial locations of zero shear on the inner and outer ring, R s3 and R s ±, respectively, in 
this case are found to be 

dV 

R s3 = 0.801 R s4 = 1.395 such that 



dz 



= where % = 3,4 (18) 

{r,z)=(R si ,0) 



We observe that the scaled poloidal shear -^== ^| z _ is negative between r = Ri and 
r = R s3 and between r R 12 and r = R s4 . Similar to the case with lids (figure 2a), 
this negative z derivative of V r means that V r < and the Ekman flow along these 
portions of rings is directed radially inward. Likewise, the positive radial velocity or 
outward Ekman flow between r = R s3 and r i? 12 and between r = R s4 and r = R 2 
corresponds to positive poloidal shear in figure 2(b). Furthermore, as in the case of lids 
(figure 2a), the signs and zeros of the scaled poloidal shear and radial velocity correlate 
well with that of the net radial force Q 2 r — |^ (dashed line in figure 2b). In addition, 
these radial locations of the reversals of net radial force and of Ekman flows near r = R s3 
and r = R sA (18) coincide within upto 2% with the radial locations R 3 and i? 4 where 
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the ideal CCF angular velocity (12) matches the angular velocity of the inner and outer 
ring Q 3 and fi 4 (table 1), namely 

R 3 = 0.793 # 4 = 1-369 such that fi< = Sl c (Ri) where i = 3,4 (19) 

This strong correlation of reversals of the net radial force with reversals of Ekman 
flow at r = R s3 and r = R sA (18), coincidental with local CCF rotation at r = R 3 w i? s3 
and r = i? 4 ps i? s4 , is completely consistent with our argument that the balance and 
unbalance between 'centrifugal' rotation and centripetal pressure gradient determines 
the fate of the radial flow along horizontal boundaries. Namely, the zero radial velocity 
at r = R s3 and r = R sA results from the CCF-like balance of centripetal pressure 
gradient — ^ and 'centrifugal' rotation Q 2 r m Vt 2 c r (12) since R 3 ~ R s3 and i? 4 m i? s4 
(18-19). Moreover, a monotonic decrease of VLq (12) with increase of r (solid line 
in figure 1) means that the near- wall fluid rotation at angular velocities of the rings fl 3 
and fi 4 is locally slower (faster) than that of CCF for the radius r that is smaller (bigger) 
than r ps R s3 and r m R s4 , correspondingly. Thus near-wall fluid rotation momentum 
deficit (excess) results, respectively, in the negative (positive) sign of the net radial force 
Q 2 r — ^ and therefore, negative (positive) sign of radial velocity V r in figure 4(b) and 
poloidal shear ^| z _ Q in figure 2(b) for the radial location r that is smaller (larger) 
than r ps R s3 and r ps _R s4 . 

In summary, the angular velocities of inner and outer rings (Q 3 and f2 4 ) set the 
CCF-like equilibrium radii (r ps R s3 and r ps i? s4 ) by matching locally to monotonically 
decreasing CCF-like profile of bulk flow rotation. The near-wall fluid over the portions 
of the rings that have a smaller radius r than these CCF equilibrium radii experience 
rotation momentum deficit that results in the inward Ekman flows due to locally higher 
centripetal pressure gradient set by faster bulk rotation as in the cases with lids. 
Conversely, when r > R s3 and r > i? s4 , the bulk rotation is slower than the near-wall 
velocity due to monotonic decrease of velocity profile with increase of radius outside 
the Ekman layers, and the fluid has enough near-wall rotation momentum to overcome 
centripetal pressure gradient and to drive the outward Ekman flows as observed in 
figure 4(b). 

The rest of the prominent features of the flow field in figure 4(b) like the strong 
vertical jets near r = R 12 and a weak outward radial jet near the midplane z = H/2 
are the direct consequences of these alternating inward-outward Ekman flows along the 
rings. Namely, driven by rotation momentum excess and deficit of fluid near inner and 
outer ring, respectively, pairs of opposing Ekman flows along both horizontal boundaries 
merge near the boundary between inner and outer ring r = R\2- Owing to continuity (6), 
these pairs of colliding Ekman flows with, presumably, equal linear radial momentum, 
launch the opposing vertical jets near the ring boundary r = R 12 that become unsteady 
with the increase of Reynolds number and break into vortex pairs or vortex rings. 
Similarly, the Ekman flows along lower and upper inner rings due to the rotation 
momentum deficit are pushed into the corners with the inner cylinder and further along 
the inner cylinder until they merge near the midplane z = H/2 to form a outward radial 
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jet as in the case with lids. But contrary to the cases with lids, the outward radial 
jet now is significantly weaker owning to the fact that the effective Reynolds number 
for these flows are smaller than in the cases with lids due to the smaller characteristic 
length scale (R s3 — Ri < L) and velocity scale (Q 3 R s3 — QiRi < U) which leads to larger 
Ekman numbers E = A q L2 = U ^r^ ^ ■ 

Finally, we would like to make two following comments. First, three-dimensional 
effects appear to be negligible at these Reynolds numbers with only noteworthy 
difference of considerably shorter vertical jets near the ring boundary r = R u as 
compared to the axisymmetric cases. Second, the angular velocities of rings control 
the angle and direction of the jet near this ring boundary r = R12. In particular, when 
rings are coupled together and rotate with the outer cylinder ('lids'), the jets become 
the inward Ekman flows along lower and upper horizontal boundary so the the angle 
with radius vector is ±7r, correspondingly. When rings are decoupled and rotate with 
the angular velocities considered above (table 1), the Ekman flows collide near the ring 
boundary r = R u and launch the opposing vertical jets, i.e. the angle is ±n/2. When 
rings are coupled to the inner cylinder, we have checked that the resulting Ekman flows 
have radially outward direction due to the excess of the near-wall angular momentum 
leading to the zero angle between the jets and radius vector in accordance with the 
mechanism described above. Moreover, this angle is expected to be sensitive to the 
details of the flow in the vicinity of the ring boundary such as presence of gaps between 
rings, three-dimensionality, etc. but it is likely to be adjusted with an appropriate choice 
of angular velocities of rings shifting the equilibrium points of local CCF balance and 
thus regulating the radial extent and radial linear momentum of the Ekman flows (cf. 
Schartman 2008). In other words, the angular velocities of rings control EC through 
the net radial momentum after the collision of Ekman flows that sets the angle at which 
the jets are launched near the ring boundary r = R 12 . 



3.1.3. Summary on Horizontal Boundary Effects The CCF-like equilibrium between 
'centrifugal' rotation and centripetal pressure gradient in cylindrical annulus is 
impossible to achieve experimentally due to the presence of the noslip horizontal 
boundaries. The rotation of these boundaries with either faster inner cylinder or slower 
outer cylinder creates the Ekman boundary layers with either angular momentum excess 
or deficit, correspondingly, and results into either outward or inward Ekman flows, 
respectively, that drive EC in the annulus. The splitting of the horizontal boundaries 
into independently rotating rings sets the CCF-like equilibrium points by matching 
locally to the CCF angular velocity, and resulting angular momentum deficit or excess 
leads to the, correspondingly, inward or outward Ekman flows along the portions of the 
rings with radius, respectively, smaller or larger than the radius of these equilibrium 
points. The opposing Ekman flows along the rings collide near ring boundaries and 
launch the vertical jets at an angle presumably determined by the mismatch of their 
linear radial momentum. This angle is expected to be sensitive to the details of the flow 
structure inside and immediately near the gaps between rings, the vertical alignment of 
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Figure 5: Dimensional azimuthal velocity profile versus radius at z = Q for experimental 

data (■) at Re ~ 9300 (Schartman 2008), circular Couette flow ( ), and 

our numerical simulations at Re = 6200: three-dimensional (— — ), axisymmetric 

( ) and noisy axisymmetric ( — • — ) cases. Three-dimensional effects are negligible 

compared to axisymmetric case both being slightly lower than experimental profile, and 
the best fit is achieved in the axisymmetric case with random noise perturbations applied 
to the surface of inner cylinder. 



the horizontal surfaces of the rings, etc. and can be adjusted by changing the angular 
velocities of the rings (cf. Schartman 2008) 

3.2. Comparison with Experiment 

We have collaborated with Princeton MRI liquid gallium experiment group and 
conducted a comparison of our computations with their experimental results. Figure 5 
shows the comparison of our numerical results for time-averaged azimuthal velocity 
in the case with lids at Re=6200 with ideal CCF profile (dotted) and experimental 
measurements (squares) conducted by Schartman (2008) at Re ~ 9300. The solid line 
corresponds to the axisymmetric computation while the fully three-dimensional results 
are shown with the dashed line. We observe that at this Reynolds number {Re ~ 6200) 
three-dimensional time-averaged azimuthal velocity is very close to the axisymmetric 
one both being upto 15% lower than the experimental data. The difference in Reynolds 
number is expected to play only a minor role in this discrepancy. 

The best fit of our (axisymmetric) computations (dash-dot line) with experimental 
data was realized when the boundary conditions (9) were perturbed with uniform 
random noise. The amplitude of the noise was 5% relative to the corresponding 
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maximums of axisymmetric solution without noise. The noise perturbation was applied 
for the part of the computational domain boundary of one spectral element long such 
that R± < r < R\+L g where L g = 0.0200 which is less than twice of the non-dimensional 
width of the gap between the inner cylinder and inner ring in the Princeton experiment 
equal to 0.0114 or 1.5 mm (Schartman 2008). 

The rational behind the noise perturbation of the boundary conditions was an 
attempt to model an effect of the centrifugally unstable flow in the gap between inner 
cylinder and inner ring. By an accident, boundary conditions on the inner cylinder 
surface (r = i?x) were also perturbed in this computation which turned out to be the 
best fit with experimental data. The effect of this perturbation of the inner cylinder 
surface boundary condition may be similar to a random blowing. This leads us to 
believe that a combination of the effects due to a run-out of inner cylinder and due 
to the centrifugally unstable flow in the gaps between the inner cylinder and inner ring 
may explain the discrepancy between the simulation and experiment (see also Schartman 
2008). 

The further work on comparison of the simulation with experiment is ongoing, and 
addition effort is needed to sort out the effects of run-out, centrifugally unstable flow 
in the gaps between cylinders and rings, vertical misalignment of horizontal surfaces of 
the rings, etc. 

3. 3. Torque and Angular Momentum Transport 

Also we have studied carefully the torque behaviour and associated angular momentum 
transport in the hydrodynamical setup of Princeton MRI liquid gallium experiment 
(Schartman 2008) as a baseline cases for our study of magneto- rotational instability 
(MRI) and MRI-driven turbulence. In order for the MRI experimental results to have 
a clear interpretation, the negative effects of the EC has to be minimized, and the 
torque amplification over the CCF torque Tq (15) with the increase of magnetic field 
can be linked directly to MRI enhancement of angular momentum transport (AMT). 
Therefore, the understanding of the torque behaviour and AMT in the baseline cases of 
hydrodynamical flow can not be overstated. 

An application of torque to the inner cylinder results in a flow that transports the 
angular momentum outward and attempts to reach a shear-free solid body rotation with 
a constant angular velocity. If the angular velocities of other boundaries are different 
form that of the inner cylinder, the resulting shear has to be maintained by application 
of torques to the boundaries in order to keep the rotation rates steady. In the context 
of MRI study, the primary interest is in the transport of the angular momentum from 
the inner to outer cylinder in centrifugally stable regime and, as shown below, in the 
minimization of EC and hence, in the reduction of the net contribution of torques exerted 
on the horizontal boundaries. Since the sum of all torques applied to the boundaries 
reflects the time increase of interior angular momentum, this contribution from the 
horizontal boundaries is equal to the sum of torques applied to the cylinders in a steady 
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Figure 6: The magnitudes of normalized torques applied to inner cylinder (Tj = T\/Tc, 
open symbols) and outer cylinder (T Q = —T 2 /Tc, filled symbols) versus Reynolds number 
for the cases with rings ( — Q— or ♦) and with lids (• -A- • • or A) in axisymmetric cases 
while the three-dimensional data are plotted with large symbols. Being the sum of 
magnitudes of torques applied to the inner and outer cylinder, the net torque exerted on 
the horizontal boundaries in the case of rings is significantly less than that in the case 
of lids making the former closer approximation to the CCF for which the net torque on 
horizontal boundaries is zero. 



state and has to be minimized for the successful MRI experiment. Note that in a case 
of unsteady flow, the time averaged torques are used when the statistically steady state 
is reached. 

The figure 6 shows the Reynolds number dependence of magnitudes of steady /time- 
averaged torque relative to the ideal CCF torque T~c (15) that has to be applied on inner 
and outer cylinders, Tj = ^ (open symbols) and T Q = (filled symbols), respectively, 
in order to keep constant boundary angular velocities (table 1). Being equal to unity for 
the case with periodic boundary conditions (CCF), the torque magnitudes are shown for 
the cases with rings (open diamonds with solid lines and filled diamonds) and lids (open 
triangles with dotted lines and filled triangles). All results are obtained in axisymmetric 
computations except for the data plotted with large open symbols that shows the results 
of fully 3D computations. Note that being zero in the ideal CCF, the difference between 
the torques applied to the inner and outer cylinder shown by open and closed symbols, 
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(a) Lids (b) Rings 

Figure 7: Steady state contour lines of effective angular momentum flux function \l> 
for the case of Re = 620 with lids in the range from -3.30 to 3.30 in the increment of 
0.31 (a) and with rings in the range from -1.35 to 1.35 in the increment of 0.21 (b). In 
the case of lids, most of the flux lines that originate from inner cylinder terminate at 
horizontal boundaries as opposed to the case of rings where they end up mostly at the 
outer cylinder which is similar to the CCF angular momentum transport between the 
cylinders. 

respectively, corresponds to the sum of torques exerted on the fluid next to the horizontal 
boundaries 

T t -T = (T 1 + T 2 )/T c 

due to zero net torque in steady/statistically steady state. Evidently, the setup 
with rings has an advantage of a smaller contribution to the net torque from the 
horizontal boundaries and of a smaller difference between inner and outer cylinder torque 
magnitudes over the setup with lids where EC is undisturbed. We also observe that in 
the range of Reynolds numbers considered the flow makes a transition from steady 
axisymmetric solution at Re = 620 to the unsteady one at Re = 6200 with small three- 
dimensional effects. Being more significant in the case with lids, three-dimensionality is 
expected to play an increasing role with the further increase of Reynolds number. 

In order to illustrate a spacial variations of AMT, we have computed an 
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effective angular momentum flux function defined in Appendix B by analogy with a 
streamfunction. The contours of the effective flux function shows the (flux) lines along 
which the angular momentum is transported, and the difference between the values of 
the flux function at two points gives the total flux across the segment of conical or 
cylindrical surfaces on which these points lie. Figure 7 shows steady state contour lines 
of constant increment for effective angular momentum flux function \l/ (B.12) for the 
case of Re = 620 with lids (a) and rings (b). For comparison, we note that the flux lines 
of (purely viscous) AMT for the ideal CCF (12) are the straight lines from inner to outer 
cylinder along z =const. Despite the fact that in both cases only a single line of 4> = 
(i.e. line of symmetry) is the same as in CCF case, the case with rings exhibits the similar 
transport of angular momentum along the flux lines that mostly originate at the inner 
cylinder and terminate at the outer cylinder in contract to the termination of the flux 
lines at the lids. The latter indicates that in the cases with lids, the angular momentum 
transport is mostly between the inner cylinder and the horizontal boundaries contrary 
to more desirable CCF-like transport between the cylinders observed in the cases with 
rings. Also note that the similarity between the shape of the flux lines away from the 
boundaries in figure 7 and the shape of vorticity contour lines and poloidal vector lines 
in figure 4 can be explained through creation of strong V r and V z components of velocity 
due to Ekman flows which affects AMT flux through advective contributions F® z and 
F% z , respectively, given by relations (B.5) and (B.6). 

In summary, if the ultimate objective is to achieve the flow with AMT as close to 
the ideal CCF as possible, the design with rings seems to have an advantage over the 
setup with lids. 

4. Conclusion and Future Work 

In this paper we have presented axisymmetric and fully three-dimensional Navier- 
Stokes calculations of circular Couette flow (CCF) in a cylindrical annulus as the first 
step in our study of magneto-rotational instability (MRI) and MRI-driven turbulence. 
Inspired by Princeton MRI liquid gallium experiment, we have computed the flow 
field in their experimental setup for realistic horizontal boundary conditions of 'lids' 
and 'rings' with the increase of Reynolds number through the onset of unsteadiness 
and three-dimensionality. The presented analysis of the flow field and angular 
momentum transport (AMT) allowed us to propose an explanation of the mechanism 
that determines the fate of the boundary flows and Ekman circulation (EC) as a result 
of a competition between the effects of 'centrifugal' rotation and pressure gradient set 
by rotation of, respectively, horizontal surfaces and bulk of the flow. In particular, with 
the appropriate choice of rotation rates of the horizontal rings that control an angle 
at which the vertical jets are launched near the ring boundaries, EC can be greatly 
diminished and CCF-like flow can be restored being more appropriate for the further 
experimental studies of MRI saturation and enhanced AMT. In addition, our numerical 
results compare favourably with the experimental data with the maximum deviation 
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below 15% being considerably smaller in the cases with 'noisy' boundary conditions. 
The future work, therefore, should involve higher Reynolds number computations with 
even more detailed modelling of experimental geometry that includes among others the 
effects of run-out of the inner cylinder, finite gaps between the cylinders and rings, and 
vertical misalignment of horizontal surfaces of the rings. 
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Appendix A. Traditional Explanation of Ekman Circulation 

Here we would like to comment that the traditional explanation of Ekman flows and 
circulation involves a balance between Coriolis and viscous forces in Ekman layers along 
rotating stressed boundary as viewed from a uniformly rotating reference frame (see 
Batchelor 1967, Greenspan 1968). As we show below, this balance of Coriolis and 
viscous forces holds, for instance, in the case of a bulk flow outside the Ekman layers 
that is close to the state of solid body rotation at, say, angular velocity of Q 2 , when the 
centripetal pressure gradient and centrifugal forces cancel each other: 

P ~lny rn 2 2 + ^~o (a.i) 

2 2 2 Or v ' 

(to the order of Ekman number E = A ^ L z ~ see Greenspan 1968) and Coriolis 
force is solely responsibly for the Ekman flow. In other flows like those considered in 
this paper, a contribution of the 'centrifugal' rotation to the Ekman flow does not reduce 
solely to Coriolis force but also includes the effects of centrifugal forces. Indeed, for the 
flow description in the non-inertial reference frame that rotates with constant angular 
velocity of the noslip boundary of the 'lids' f2 2 = ^2 e z , the centrifugal and Coriolis 
body forces has to be added to the right-hand side of equations (3-5), 

f — — Q, 2 x (f2 2 x t) — 2 f2 2 x u 

f r = tt 2 2 r + 2 tt 2 u e (A.2) 

along with the replacement of the inertial frame velocity V with the velocity in the non- 
inertial rotating reference frame u leaing to the following radial momentum equation in, 
e.g., case of axisymmetry (^ = 0) 



du r „ u%. ^ 9 „ ^ 1 T a u r 

+ u • V - = ttjr + 2 n 2 u e + — A (r:Z) u r - — 
at r Re L r l 

or 



^(r,z)Uf 



dp 
dr 



9U * S1 rr + (2n 2U . e + ^-^ + -L^ + ... (A.3) 



dt \ r J dr Re dr 2 

Let us now compare this relation with equation (16) where we expand the 'centrifugal' 
term Q 2 r in powers of azimuthal velocity deviation from the state of solid body rotation 
of the lids, U = V - rVt 2 : 



Therefore, equation of radial momentum balance (16) rewritten with the expansion (A.4) 
coincides exactly with equation (A.3) in the case of axisymmetric flow when the 
difference between V r and u r and between Ue and uq disappears due to the relationship 
between velocities in the inertial and non-inertial reference frames: 

V(r, 9, z,t) = r Vt 2 e® + u(r, = 9 + Q 2 t, z, t) 



or 



V(r,z,t) =rQ 2 e +u(r,z,t) (A.5) 
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In a particular case of flows that are close to solid body rotation when considered in 
this description (A. 3) equivalent to our earlier framework (16), the contribution of the 
centrifugal and pressure gradient forces cancels out (A.l) resulting in the steady state 
balance of Coriolis and viscous forces, 

1 d 2 u 

= 2Q 2 u e + -^ + --- (A.6) 

u 2 

where we omitted higher order terms including C 2 fi 2 «e away from the axis of 
rotation r = due to vanishing Me at the noslip boundary. In more general case, the 
additional contributions of the centrifugal forces Q?,r and — to the Ekman flows has to 
be considered being a part of the 'centrifugal' term Vt 2 r (16, A. 3) that is balanced by 
centripetal pressure gradient and viscous forces in the Ekman flows along the lids that 
drive EC in the cylindrical annulus. 



Appendix B. Angular Momentum Flux and Flux Function 

Let us take the axisymmetric version of azimuthal momentum equation (4) with viscosity 
v instead of and rewrite it in terms of conservation of axial angular momentum 
C = rVg. The summation of the axisymmetric versions of equations (4) and (6) 
multiplied by factors — r and rVg, correspondingly, gives 



d_ 
dt 



[rVa 



— r 



T . dV e dV e V r V e 

-V r — V z — h V 



dr 



dz 



d 2 V e IdVp d 2 V e Vg 
dr 2 r dr dz 2 r 2 



+ rV e 



dV r | dV z | V r 
dr dz r 



rV r 



dr 



+ rV 2 



dVo_ 

dz 



+ V r V e - v { r 



+ 



dV r dV z 
rV d —^+rV e -^ + V r V e 
dr dz 
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(B.l) 
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The above equation (B.l) reflects the conservation of axial or ^-component of 
angular momentum and can be condensed to 
f) 1 f) f)F 

S^+W^ + lif =° ( R2 » 

or 

^(rx V) + V (r ,, r F) =0 (B.3) 

where total angular momentum flux tensor F and its components F rz and F zz are given 
for axisymmetric flows by 

F = F a + F v (B.4) 

F r a z = rV r Ve (B.5) 

K z = rV e V z (B.6) 

F: z = - ,r 2 |- (Y?) = -rr r9 (B.7) 

K z = - vr-£ = -rr ze (B.8) 

where superscript a and v denote an advective and viscous contributions to the total 
flux of angular momentum, and r r9 and r zd are components of the stress tensor r. 

By analogy with poloidal streamfunction ip that satisfies axisymmetric version of 
continuity equation (6) due to, e.g., a definition 

M _1 fl£« 

oz r or 

an angular momentum flux function \l/(r, z) for steady flows can be introduced to satisfy 
the steady version of equation (B.2): 

F „ = » F „ = -im. (b.io) 

oz r or 

As the lines of constant value of poloidal streamfunction ip are lines of constant poloidal 
flow rate with poloidal velocity vector being tangent to these streamlines, the lines of 
constant total angular momentum flux function ^ are the flux lines along which angular 
momentum is transported, and the difference between these values at any two points 
gives the total flux of angular momentum transferred across the line that joints these 
points. Note that in the case of unsteady flow, the time-averaged quantities can be used 
in statistically steady state instead of instantaneous ones. Also note that similar to the 
streamfunction, the flux function \1> as a solution of equations (B.10) is defined up to a 
constant. We have used the point of half height on inner cylinder as a zero value for the 
flux function 

* (r = R u z = ^ J = (B.ll) 



and for convenience, we define effective flux function f as a flux function \l> multiplied 
by circumference 2irr and scaled by the ideal CCF torque T c (15) 
~ 2ttt 

m(r,z) = —m(r,z) (B.12) 
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Finally, due to multiplication of flux function \l/ by r, the difference in effective flux 
function ^ between two points gives the total flux of angular momentum across the 
conical or cylindrical surface on which these points lie. Thus, in the case of e.g. inner 
and outer cylinder, it can be shown that torques T x and T 2 scaled by the ideal CCF 
torque Tc are given by 

J- = H) - 0) Z Z 1 = *(i*2, H) - *(i? 2 , 0) (B.13) 

Similar expressions hold for the torques applied to the horizontal surfaces. 
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1. Introduction 

The phenomenon of Ekman circulation (EC) occurs in most if not all rotating flows 
with stressed boundaries that are not parallel to the axis of rotation. The manifestation 
of EC ranges from wind-driven ocean currents (?), to the accumulation of the tea 
leaves at the bottom of a stirred cup (see, e.g., ?). One of consequences of EC and 
of the associated Ekman flows is greatly to to enhance mixing and transport and in 
particular, the transport of angular momentum, above the values due to viscosity alone. 
Traditionally, Ekman flows are explained in terms of action of Coriolis forces in the 
Ekman layers along the rotating stressed boundaries (?). 

There are circumstances when the presence of EC has undesirable effects. For 
example, this is the case in laboratory experiments to study the development of magneto- 
rotational instability (MRI) in liquid metals (see a monograph edited by ?). The MRI 
instability is important in astrophysics where it is believed to lead to turbulence in 
magnetized accretion disks (?). Many of the features of the MRI and its associated 
enhancement of angular momentum transport (AMT) can be studied experimentally in 
magnetized flows between rotating coaxial cylinders. In these experiments, the rotation 
rates of the cylinders are chosen in such a way that the fluid's angular momentum 
increases outwards so that the resulting rotational profile is stable to axisymmetric 
perturbations (so-called centrifugally stable regime). The presence of a weak magnetic 
field can destabilize the basic flow, provided the angular velocity increases inward, and 
lead to an enhancement of outward AMT. 

In an ideal situation, the basic state consists of circular Couette flow (CCF), and 
the outward transport of angular momentum in the absence of magnetic fields is solely 
due to viscous effects. The presence of a magnetic field would destabilize the basic 
flow through the effects of MRI and lead to a measurable increase of AMT. In practice, 
this ideal case can never be realized in laboratory experiments because of horizontal 
boundaries. The presence of these boundaries drives an EC that enhances AMT even 
in the absence of magnetic effects. In order to study the enhancement of AMT due to 
MRI it is crucial to be able to distinguish the effects that are magnetic in origin from 
those that are due to the EC. One possibility is to make the cylinders very tall so the 
horizontal boundaries are far removed from the central region. This approach, however, 
is not practical owing to the high price of liquid metals. 

The alternative approach is to device boundaries in such a way that the resulting 
EC can be controlled and possibly reduced. For example, attaching the horizontal 
boundaries to the inner or outer cylinder results in dramatically different flow patterns. 
Another possibility could be to have the horizontal boundaries rotating independently 
of inner and outer cylinder. Goodman, Ji and coworkers (?, ?, ?) have proposed to 
split the horizontal boundaries into two independently rotating rings whose rotational 
speeds are chosen so as to minimize the disruption to the basic CCF by secondary Ekman 
circulations. Indeed this approach has been implemented in the Princeton's MRI liquid 
gallium experiment (?). In any case, no matter how the horizontal boundaries are 
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implemented it is important to understand what kind of EC patterns arise before the 
magnetic effects are introduced. 

In the present paper we address this issue by studying the effects of horizontal 
boundary conditions on CCF numerically. We study both axisymmetric and fully three- 
dimensional geometries and investigate the effects of changing rotation rates (Reynolds 
number) through the onset of unsteadiness and three-dimensionality. The next section 
(section 2) describes the formulation of the problem and gives an account of numerical 
aspects of its solution technique including a brief description of the spectral element 
code Nek5000 (?). The section 3 starts with an explanation of flow behaviour due to 
horizontal boundary conditions, i.e. CCF, Ekman and disrupted Ekman circulation due 
to periodic horizontal boundaries, 'lids' and 'rings', correspondingly (section 3.1). Then 
the paper proceeds with description of comparison of our results with the experimental 
data (section 3.2) followed by an examination of torque and AMT (section 3.3). Finally, 
we draw conclusions and describe future work in section 4. 



2. Problem Formulation and Numerical Method 



2.1. Formulation 



We study the flow of an incompressible fluid with finite (constant) kinematic viscosity 
v in a cylindrical annulus bounded by coaxial cylinders. The cylinders have the radii 
R\ and R 2 (-^i < -^2) an d rotate with angular velocities Q* and f^, respectively. The 
annulus is confined in the vertical direction by horizontal boundaries at distance H* 
apart. The formulation of the problem in cylindrical coordinates (r, 9, z) with the scales 
for characteristic length L and velocity U, 

L = R* 2 -R\ U = QlRl - tt* 2 R* 2 (1) 

and therefore, with the relationship between dimensional variables (with asterisk) and 
non-dimensional radius, height, velocity vector V, time and pressure given by 



[r*,z*, V*,t*,p*] 



Lr,Lz,UV,-t,pU 2 p 



(2) 



correspondingly, results in the following non-dimensional incompressible Navier-Stokes 
equations: 



8V r 

~di 

m 

dt 
dt 

dV r 



+ (V- V 



V 2 1 
r Re 



V r V e 



+ (V • V) V e + — - = — 
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dr 
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I dp 
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~ rd9 



dz r 
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Re dz 
= 



dr r d9 

where p is a constant fluid density and Reynolds number Re is defined as 
UL _ (QlRl - n* 2 R* 2 )(R* 2 - R{) 
v v 



Re 



(3) 

(4) 

(5) 
(6) 

(7) 
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while the scalar advection operator due to a vector field V and laplacian of a scalar 
function S(r, z) are given by 



dS V e dS 



dS 



(V -V)S = V r — + — — + V z — AS 



r dr ' r d6 2 dz dr 2 r dr r 2 d6 2 ' dz 2 



d 2 S ldS_ 

dr 2 r dr 



1 d 2 S d 2 S 



(8) 



The initial conditions for the flow in the annulus and boundary conditions at the 
cylinder surfaces r — R\ and r = R 2 are 

V r = V z = V e = rtt(r) (9) 

where non-dimensional angular velocity Q(r) is given by circular Couette flow (CCF) 
profile 



n c (r) = A + 



B 



Q 2 R 2 — Q\R\ 



m 



B 



R^R^i^l — 1^2; 



(10) 



^3 

n 2 



R 2 — R\ 

At the horizontal boundaries z = and z = H, two types of the boundary conditions 
have been considered, namely, lids and rings, given by (9) where angular velocity Q(r) 
is equal to 

r = Ri 

R 1 <r<R l2 (n) 
R 12 <r<R 2 [ ' 

r = R 2 

Here R\ 2 is the radial location of the boundary between the inner and outer rings, and 
Q 3 and are angular velocities of inner and outer rings, correspondingly. Inspired by 
Princeton MRI liquid gallium experiment (?), the non-dimensional angular velocities 
and cylinder height as well as cylinder and ring boundary radii used in this study are 
given in table 1 in addition to the dimensional parameters involved in comparison with 
the experiment (subsection 3.2). In the cases with lids, angular velocities f2 3 and f2 4 are 
equal to the angular velocity of the outer cylinder fl 2 while in the cases with rings they 
turn out to be close to the values of CCF profile (10) taken at the middle of radii of the 
corresponding rings. 



2.2. Numerical Technique 

The axisymmetric version of equations (3-6) and fully three-dimensional version in 
cartesian coordinates has been solved numerically with the spectral-element code 
Nek5000 developed and supported by Paul Fischer and collaborators (see ?, ?, and 
references within). 

The temporal discretization in Nek5000 is based on a semi-implicit formulation 
in which the nonlinear terms are treated explicitly in time and all remaining linear 
terms are treated implicitly. In particular, we used either a combination of fcth- 
order backward difference formula (BDFA;) for the diffusive/so lenoidal terms with 
extrapolation (EXTA; — 1) for the nonlinear terms or the operator-integration factor 
scheme (OIFS) method where BDF/c is applied to the material derivative with the 
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Lids Rings 


Experiment: Lids 


Ri 


0.538 


R\ 


(cm) 


7.1 


R2 


1.538 




(cm) 


20.3 


R12 


1.038 


r?* 

K 12 


(cm) 


13.7 


H 


2.114 


H* 


(cm) 


27.9 


fix 


3.003 




(rpm) 


200 


n 2 


0.488 0.400 


n* 2 


(rpm) 


26 


^3 


0.488 1.367 


n* 3 


(rpm) 


26 


fi 4 


0.488 


n* 4 


(rpm) 


26 




Table 1: The geometry and rotation parameters for the computational cases with lids 
and rings at Re = 6190 and experimental setup with lids at Re = 9270 along with the 
drawing of the cut of 3D computational mesh at 9 = for the case with rings. Note 
the clustering of the gridlines at the boundaries of the spectral elements whose location 
and dimensions are chosen to resolve efficiently boundary layers and 'step' changes in 
angular velocity between cylinders and rings. 



explicit fourth-order Runge-Kutta scheme being used for the resulting pure advection 
initial value problem. 

With either the BDF/c/EXTA; — 1 or OIFS approach, the remaining linear portion 
of time advancement amounts to solving an unsteady Stokes problem. This problem 
is first discretized spatially using spectral-element method (SEM) and then split into 
independent subproblems for the velocity and pressure in weak variational form. The 
computational domain is decomposed into K non-overlapping subdomains or elements, 
and within each element, unknown velocity and pressure are represented as the tensor- 
product cardinal Lagrange polynomials of the order N and N — 2, correspondingly, 
based at the Gauss-Lobatto-Legendre (GLL) and Gauss-Legendre (GL) points. This 
velocity-pressure splitting and GLL-GL grid discretization requires boundary condition 
only for velocity field and avoids an ambiguity with the pressure boundary conditions 
in accordance with continuous problem statement. 

The discretized Stokes problem for the velocity update gives a linear system which 
is a discrete Helmholtz operator. It comprises the diagonal spectral element mass matrix 
with spectral element Laplacian being strongly diagonally dominant for small timesteps, 
and therefore, Jacobi (diagonally) preconditioned conjugate gradient iteration is readily 
employed. Then the projection of the resulting trial viscous update on divergence-free 
solution space enforces the incompressibility constraint as the discrete pressure Poisson 
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equation is solved by conjugate gradient iteration preconditioned by either the two- level 
additive Schwarz method or hybrid Schwarz/multigrid methods. Note that we used 
dealising/overintegration where the oversampling of polynomial order by a factor of 
3/2 was made for the exact evaluation of quadrature of inner products for non-linear 
(advective) terms. 

The typical axisymmetric case with rings at high Reynolds number of Re = 6200 
(see figure 4b) required the spacial resolution with polynomial order N — 10 and number 
of spectral elements K = 320 (cf. drawing for table 1) and was computed with timestep 
At = 10~ 3 for the duration of t ~ 300, while the axisymmetric run with lids at the same 
Re (figure 4a) had N = 8, K = 476, At = 5 x 10~ 3 and t ~ 500. The corresponding 
three-dimensional cases with rings and lids had N — 11, K — 9600, At = 6.25 x 10~ 4 , 
t ~ 280 and N = 9, K = 14280, At = 6.25 x 10~ 4 , t ~ 180, respectively. Note that in 
order to facilitate time advancement and minimize CPU requirements, the final output 
from another cases, e.g. with lower Reynolds number Re, was used as initial conditions 
for some of the computations with higher Re, and the corresponding axisymmetric 
cases with small random non-axisymmetric perturbation was a starting point for most 
of our fully 3D computations. Apart from CPU savings, the usage of the perturbed 
axisymmetric solution obtained in cylindrical formulation (3-6) as initial condition for 
3D computations at low Reynolds numbers (Re = 620) served as an additional validation 
of the code setup due to the convergence of the fully 3D results computed in cartesian 
formulation back to the unperturbed axisymmetric steady state initial condition (see 
also subsection 3.3). 

Finally, the step change of angular velocities that mimics its transition in the gaps 
or grooves between the cylinders and horizontal boundaries as well as between the inner 
and outer ring in Princeton MRI liquid gallium experiment (?) was modelled within one 
spectral element of the radial size L g = 0.020 by ramping power law function of radius 
with an exponent that was varied in the range from 4 to iV — 1 without noticeable effect 
on the flow. 

3. Results 

Let us first start with examination of effects of horizontal boundary conditions on flow 
pattern in general and Ekman circulation in particular before moving to a comparison 
with the experiment and examination of angular momentum transport in the cylindrical 
annulus. 

3.1. Horizontal Boundary Effects 

Here we contrast two type of horizontal boundary conditions with an ideal baseline case 
of circular Couette flow (CCF). Being zero in the ideal case, we argue that the unbalance 
between 'centrifugal' rotation and centripetal pressure gradient determines the fate of 
the radial flow along horizontal boundaries in the cylindrical annulus. 
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r 

Figure 1: The CCF azimuthal velocity Vg (- -), angular velocity Qc {~ ~~ ) ; axial 

angular momentum C = rVg ( — ■ — ) and minus pressure Pq ( ) versus radius. 

Note the monotonically increasing angular momentum and decreasing angular velocity 
with radius for centrifugally stable circular Couette flow where the 'centrifugal' rotation 
balances the centripetal pressure gradient leading to zero radial and axial velocities. 



In the ideal case of CCF, the sheared circular motion is balanced by centripetal 
pressure gradient. To be precise, the ideal CCF is the following exact solution of 
equations (3-6) for periodic (or stress-free) horizontal boundary conditions: 

V r = V z = V e = rQ c (r) = Ar + - 

r 

rr y2 £2 2 g2 

Pc(r) = / — dr = -— - — + 2AB log r + Const (12) 
J r 2 2r z 

Here the constant A given by equation (10) is proportional to the increase in axial 
angular momentum, 

C = rVg = Q r 2 (13) 

outward between the cylinders while the constant B is set by shear-generating angular 
velocity drop between them. The figure 1 shows CCF azimuthal velocity V e (dashed), 
angular velocity Qc (solid), axial angular momentum £ (dash-dotted) and negative of 
pressure, —pc (dotted,) for the non-dimensional parameters given in table 1. Since we 
are primarily interested in further MRI studies, the baseline flow has to be centrifugally 
stable, i.e. with angular momentum C increasing outward (for Q > 0), and therefore, 
satisfying Rayleigh criterion 

dC 2 

4r > ° (") 

which is the case in this ideal CCF (dash-dotted line in figure 1). In order to maintain 
rotation with shear in this virtual experiment with periodic horizontal boundaries, the 
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positive axial torque Tq 

-H f 2iT a 



r = /( ? x(ri.r)) f =/ / 



r d V$ 

dzdO— 

He or r 



A7rHR 2 R 2 Q 1 -Q 2 

(15) 



V e =rn c "^2 — -^1 R e 

has to be applied to the inner cylinder while the outer cylinder is kept from shear-free 
solid body rotation (0,(r) = fii = A, B = T = 0) by negative torque —T c . Note that in 
the above equation (15), r is the non-dimensional shear stress tensor (see also Appendix 
B). 

3.1.1. Ekman Circulation with 'Lids' In practice, the ideal CCF can never be realized 
in laboratory experiments because of horizontal boundaries. The simplest realizable 
configuration is the one we refer to as 'lids' when horizontal boundaries are coupled to 
the outer cylinder (fi 3 — Q 4 — Q 2 )- To see how flow changes in the presence of lids that 
rotate with outer cylinder, let us imagine that these lids were inserted impulsively into 
fluid with ideal CCF profile given by equation (12) and plotted as a solid line in figure 1 
for Qi and Q 2 from table 1. Keeping the most important terms in the axisymmetric 
form of the equation (3) gives 

dV r Q2 dp 1 d 2 V r 

-dT = nr -d-r + R-e^ + --- (16) 
where we used Vg = rfl. For the initial condition of CCF (12), the left-hand side of 
equation (16) is equal to zero everywhere outside the lids which is also consistent with 
zero radial flow V r — 0. This zero radial flow also results in zero diffusion term -^^pf in 
equation (16) and zero net radial force Vl 2 r — The latter results from the exact CCF 
balance between (positive) 'centrifugal' rotation term fl^ r and (negative) centripetal 
pressure gradient term — |^ in equation (16). 

Instead of initial ideal CCF angular velocity Qc (12), the flow next to the lids now 
rotates with a smaller angular velocity of the outer cylinder (Q 2 = ^3 = ^4 < ^c)- 
However, the centripetal pressure gradient is still set by the bulk rotation of the rest of 
the fluid and therefore, becomes suddenly larger than the 'centrifugal' rotation of fluid 
next to the lids, i.e. |^ = Vt^r > VL 2 r. As a result of this angular momentum deficit 
of near- wall fluid, the centripetal pressure gradient prevails over rotation term in (16). 
Therefore, the net radial force becomes non-zero and negative, Vt 2 r — |^ < 0, resulting 
in negative sign of ^ (16) and therefore, in a formation of the Ekman layer with an 
inward radial flow (V r < 0) in the vicinity of the lids. 

The figure 2(a) confirms that the net radial force near, e.g. the lower horizontal 
surface z — 0, Q 2 r — §f | 2 _ (dashed) is negative, as well as the scaled poloidal wall 
shear -2= ^r-M „. The latter means that the z-derivative of V r is negative at the lower 

VRe oz 1 2=0 7 & 

lid which in turn results in a decrease of radial velocity with the increase of height z 
from noslip zero value at the lid, V r = (9) to negative values associated with the 

z=0 

inward Ekman flow. Thus the deficit of angular momentum in the near-wall fluid of the 
Ekman layer results in unbalanced centripetal pressure gradient set by the bulk rotation 
of the rest of the flow outside the layer and drives the Ekman flow radially inward. 
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(b) Rings 

Figure 2: Steady state scaled radial wall shear (— — ) and near- wall net radial force 
(- -) for Re = 620 in the case of lids (a) and rings (b). The definite negative net 
radial force in the lids case (a) result in the inward radial Ekman flow with negative 
radial wall shear being disrupted in the case of rings (b) by the alternating sign of the 
net radial force that correlates well with the sign of the wall shear and thus with the 
alternating directions of Ekman flows. 
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To summarize, the presence of slower rotating lids disrupt the initial ideal CCF 
equilibrium between centrifugal rotation and centripetal pressure gradient set by 
rotation of, respectively, lids and bulk of the flow. This leads to the negative net radial 
force, Vt 2 r — ^ < and inward Ekman flow, V r < owning to ^ < 0. As time grows, so 
does the magnitude of negative radial velocity in the Ekman layer and, eventually, the 
diffusion term ^^Jr (16) in the Ekman boundary layer of the width Az ~ 0{\/R~e) 
becomes of the same order (i.e. ~0(1)) as the net radial force that results from two 
other terms in (16). Thus the diffusion effects finally balance the rotation momentum 
deficit of the fluid in Ekman boundary layers near the lids in the saturation steady state 
(see also Appendix A). 

To check consistency of this argument, the saturation magnitude of || across the 
layer is verified to be more than an order of magnitude smaller than the corresponding 
1^ which confirms that saturation centripetal pressure gradient — 12 is indeed set by 
the bulk rotation of the fluid outside the Ekman boundary layers. The saturation 
bulk rotation can be illustrated by the instantaneous saturation profiles of azimuthal 
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Figure 3: Azimuthal velocity Vq versus radius for the circular Couette flow (— — ) and 
instantaneous axisymmetric profiles at z = -j in the case of lids for the series of Reynolds 

numbers Re = 620 ( ), 1900 (dash-triple-dot) and 6200 ( — • • — ), and in the case 

of rings for the same Reynolds numbers: ( ), ( — • — ) and (- -), respectively. 

The Ekman- circulation induced momentum deficiency in azimuthal velocity profiles in 
the cases with lids is greatly diminished by the particular choice of angular velocities of 
independently rotating rings. 



velocity shown in figure 3 for the steady cases with lids for Re = 620 (dashed) and 
1800 (dash-triple-dot), and unsteady case with lids of Re = 6200 (dash-double-dot) 
at z = ^j. Interesting that for the range of Reynolds numbers considered, the effect 
of the increase of Reynolds number is minor in comparison with significant azimuthal 
momentum deficiency resulted from the change of horizontal boundary conditions from 
initial ideal CCF (solid) to the cases of Ekman flows over lids. 

Owing to momentum deficiency of the near wall fluid, the higher centripetal 
pressure gradient drives the inward Ekman flows along the lids that result in EC in 
the cylindrical annulus. In order to further illustrate the phenomenon of EC due to 
horizontal boundaries, we have plotted the contours of azimuthal vorticity ujq and vector 
plot of poloidal velocity (V r) V z ) in figure 4(a) for the case of Re = 6200 with the former 
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(a) Lids 




(b) Rings 



Figure 4: Instantaneous contours of azimuthal vorticity and vector field of poloidal 
velocity for Re = 6200 in the case of lids (a) and rings (b). The Ekman circulation and 
outward radial jet near the midplane in the case of lids is severly disrupted in the setup 
with rings due to alternating inward- outward Ekman flows. 



given by 



dV r _ dV z 
dz dr 



Ug 



8V r 

dz 



LUg 



~-Ri 



dr 



(17) 



where noslip conditions (9) along the walls has been used. Here the vorticity contours 
are coloured from blue {uj e < 0) to red {ujg > 0). Note that this change of colours from 
blue to red (through green colour whenever it is visible) shows the locus of zero vorticity 
that gives approximate location of a jet or jet-like features in the flows along the lids at 
z = or z = H with extremum in V r and in the flows along the cylinders at r = R% or 
r = R 2 with minimum or maximum in V z (17). In figure 4(a), we observe the Ekman 
boundary layers along the lids with vorticity contours changing their colours from blue 
{ujq < 0) to red at the lower lid and from red (ug > 0) to blue at the upper lid. In both 
instances, this change of colour shows the locus of minimum in V r < or location where 
inward Ekman flow is the strongest. Similarly, the change of colours near the inner 
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cylinder surface shows the opposing vertical jet-like flows along the inner cylinder that 
merge near the midplane z = H/2 and owing to continuity (6), form strong outward 
radial jet. At these high Reynolds numbers, beyond Re ~ 1800, the radial jet becomes 
unsteady and starts to oscillates breaking into pairs of vortices or to be precise, into 
pairs of vortex rings that move toward the lids drawn by the mass loss in the Ekman 
layers and thus closing the EC cycle. 

Summing up the flow pattern in the case of lids, we conclude that because of the 
deficit of rotation momentum in Ekman layers, the fluid is pushed centripetally in these 
layers along the lids and further along the inner cylinder with the subsequent formation 
of the strong outward radial jet that eventually transport fluid back to the lids and 
closes the cycle of the EC (see also Appendix A and section 3.3). 

3.1.2. Ekman Circulation Disruption Due to 'Rings' When each horizontal boundary 
is split into a pair of rings that rotate independently with the angular velocities Q 3 and 
Q 4 (table 1), the bulk rotation and resulting centripetal pressure gradient is restored 
back to that of the CCF. The restored profiles of azimuthal velocity in the cases with 
rings are shown in figure 3 for the same Reynolds numbers as for the cases with lids, 
namely for the steady cases of Re = 620 (dotted) and 1800 (dash-dot) and unsteady 
case of Re = 6200 (long dash). Along with the restoration of the bulk rotation back to 
that of the CCF, we observe other major differences between the cases with lids (a) and 
rings (b) in the flow field structure illustrated in figure 4. 

Instead of a single outward radial jet and inward Ekman flows along the lids, 
figure 4(b) shows alternating inward-outward Ekman flows along the rings that, as we 
describe below, produce strong vertical jets near r = R\ 2 and a weaker outward radial 
jet near midplane z = H/2. The alternating inward-outward Ekman flows along, e.g., 
the lower inner and outer rings [z — 0), are also evident in figure 2(b) where the scaled 
poloidal shear ^= ^| z=0 (solid) is plotted as a function of radius for Re = 620. The 
radial locations of zero shear on the inner and outer ring, R s3 and R s4 , respectively, in 
this case are found to be 

8Vr 

R s3 = 0.801 R sA = 1.395 such that 



dz 



= where % = 3,4 (18) 

{r,z)=(R 3i ,0) 



We observe that the scaled poloidal shear ^ | ___ Q is negative between r — Ri and 
r = R s3 and between r m R12 and r = R s4 . Similar to the case with lids (figure 2a), 
this negative z derivative of V r means that V r < and the Ekman flow along these 
portions of rings is directed radially inward. Likewise, the positive radial velocity or 
outward Ekman flow between r = R s3 and r Ru and between r = R s4 and r = R 2 
corresponds to positive poloidal shear in figure 2(b). Furthermore, as in the case of lids 
(figure 2a), the signs and zeros of the scaled poloidal shear and radial velocity correlate 
well with that of the net radial force Q 2 r — |^ (dashed line in figure 2b). In addition, 
these radial locations of the reversals of net radial force and of Ekman flows near r = R s3 
and r = R s4 (18) coincide within upto 2% with the radial locations R 3 and R 4 where 
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the ideal CCF angular velocity (12) matches the angular velocity of the inner and outer 
ring Q 3 and fi 4 (table 1), namely 

R 3 = 0.793 # 4 = 1-369 such that fi< = Sl c (Ri) where i = 3,4 (19) 

This strong correlation of reversals of the net radial force with reversals of Ekman 
flow at r = R s3 and r = R sA (18), coincidental with local CCF rotation at r = R 3 w i? s3 
and r = i? 4 ps i? s4 , is completely consistent with our argument that the balance and 
unbalance between 'centrifugal' rotation and centripetal pressure gradient determines 
the fate of the radial flow along horizontal boundaries. Namely, the zero radial velocity 
at r = R s3 and r = R sA results from the CCF-like balance of centripetal pressure 
gradient — ^ and 'centrifugal' rotation Q 2 r m Vt 2 c r (12) since R 3 ~ R s3 and i? 4 m i? s4 
(18-19). Moreover, a monotonic decrease of VLq (12) with increase of r (solid line 
in figure 1) means that the near- wall fluid rotation at angular velocities of the rings fl 3 
and fi 4 is locally slower (faster) than that of CCF for the radius r that is smaller (bigger) 
than r ps R s3 and r m R s4 , correspondingly. Thus near-wall fluid rotation momentum 
deficit (excess) results, respectively, in the negative (positive) sign of the net radial force 
Q 2 r — ^ and therefore, negative (positive) sign of radial velocity V r in figure 4(b) and 
poloidal shear ^| z _ Q in figure 2(b) for the radial location r that is smaller (larger) 
than r ps R s3 and r ps _R s4 . 

In summary, the angular velocities of inner and outer rings (Q 3 and f2 4 ) set the 
CCF-like equilibrium radii (r ps R s3 and r ps i? s4 ) by matching locally to monotonically 
decreasing CCF-like profile of bulk flow rotation. The near-wall fluid over the portions 
of the rings that have a smaller radius r than these CCF equilibrium radii experience 
rotation momentum deficit that results in the inward Ekman flows due to locally higher 
centripetal pressure gradient set by faster bulk rotation as in the cases with lids. 
Conversely, when r > R s3 and r > i? s4 , the bulk rotation is slower than the near-wall 
velocity due to monotonic decrease of velocity profile with increase of radius outside 
the Ekman layers, and the fluid has enough near-wall rotation momentum to overcome 
centripetal pressure gradient and to drive the outward Ekman flows as observed in 
figure 4(b). 

The rest of the prominent features of the flow field in figure 4(b) like the strong 
vertical jets near r = R 12 and a weak outward radial jet near the midplane z = H/2 
are the direct consequences of these alternating inward-outward Ekman flows along the 
rings. Namely, driven by rotation momentum excess and deficit of fluid near inner and 
outer ring, respectively, pairs of opposing Ekman flows along both horizontal boundaries 
merge near the boundary between inner and outer ring r = R\2- Owing to continuity (6), 
these pairs of colliding Ekman flows with, presumably, equal linear radial momentum, 
launch the opposing vertical jets near the ring boundary r = R 12 that become unsteady 
with the increase of Reynolds number and break into vortex pairs or vortex rings. 
Similarly, the Ekman flows along lower and upper inner rings due to the rotation 
momentum deficit are pushed into the corners with the inner cylinder and further along 
the inner cylinder until they merge near the midplane z = H/2 to form a outward radial 
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jet as in the case with lids. But contrary to the cases with lids, the outward radial 
jet now is significantly weaker owning to the fact that the effective Reynolds number 
for these flows are smaller than in the cases with lids due to the smaller characteristic 
length scale (R s3 — Ri < L) and velocity scale (Q 3 R s3 — QiRi < U) which leads to larger 
Ekman numbers E = A q L2 = U ^r^ ^ ■ 

Finally, we would like to make two following comments. First, three-dimensional 
effects appear to be negligible at these Reynolds numbers with only noteworthy 
difference of considerably shorter vertical jets near the ring boundary r = R u as 
compared to the axisymmetric cases. Second, the angular velocities of rings control 
the angle and direction of the jet near this ring boundary r = R12. In particular, when 
rings are coupled together and rotate with the outer cylinder ('lids'), the jets become 
the inward Ekman flows along lower and upper horizontal boundary so the the angle 
with radius vector is ±7r, correspondingly. When rings are decoupled and rotate with 
the angular velocities considered above (table 1), the Ekman flows collide near the ring 
boundary r = R u and launch the opposing vertical jets, i.e. the angle is ±n/2. When 
rings are coupled to the inner cylinder, we have checked that the resulting Ekman flows 
have radially outward direction due to the excess of the near-wall angular momentum 
leading to the zero angle between the jets and radius vector in accordance with the 
mechanism described above. Moreover, this angle is expected to be sensitive to the 
details of the flow in the vicinity of the ring boundary such as presence of gaps between 
rings, three-dimensionality, etc. but it is likely to be adjusted with an appropriate choice 
of angular velocities of rings shifting the equilibrium points of local CCF balance and 
thus regulating the radial extent and radial linear momentum of the Ekman flows (cf. 
?). In other words, the angular velocities of rings control EC through the net radial 
momentum after the collision of Ekman flows that sets the angle at which the jets are 
launched near the ring boundary r = i? 12 . 



3.1.3. Summary on Horizontal Boundary Effects The CCF-like equilibrium between 
'centrifugal' rotation and centripetal pressure gradient in cylindrical annulus is 
impossible to achieve experimentally due to the presence of the noslip horizontal 
boundaries. The rotation of these boundaries with either faster inner cylinder or slower 
outer cylinder creates the Ekman boundary layers with either angular momentum excess 
or deficit, correspondingly, and results into either outward or inward Ekman flows, 
respectively, that drive EC in the annulus. The splitting of the horizontal boundaries 
into independently rotating rings sets the CCF-like equilibrium points by matching 
locally to the CCF angular velocity, and resulting angular momentum deficit or excess 
leads to the, correspondingly, inward or outward Ekman flows along the portions of the 
rings with radius, respectively, smaller or larger than the radius of these equilibrium 
points. The opposing Ekman flows along the rings collide near ring boundaries and 
launch the vertical jets at an angle presumably determined by the mismatch of their 
linear radial momentum. This angle is expected to be sensitive to the details of the flow 
structure inside and immediately near the gaps between rings, the vertical alignment of 
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Figure 5: Dimensional azimuthal velocity profile versus radius at z = x for experimental 

data (■) at Re ~ 9300 (?), circular Couette flow ( ), and our numerical simulations 

at Re = 6200: three-dimensional (— — ), axisymmetric ( ) and noisy axisymmetric 

( — • — ) cases. Three-dimensional effects are negligible compared to axisymmetric case 
both being slightly lower than experimental profile, and the best fit is achieved in the 
axisymmetric case with random noise perturbations applied to the surface of inner 
cylinder. 



the horizontal surfaces of the rings, etc. and can be adjusted by changing the angular 
velocities of the rings (cf. ?) 

3.2. Comparison with Experiment 

We have collaborated with Princeton MRI liquid gallium experiment group and 
conducted a comparison of our computations with their experimental results. Figure 5 
shows the comparison of our numerical results for time-averaged azimuthal velocity 
in the case with lids at Re=6200 with ideal CCF profile (dotted) and experimental 
measurements (squares) conducted by Schartman (2008) at Re ~ 9300. The solid line 
corresponds to the axisymmetric computation while the fully three-dimensional results 
are shown with the dashed line. We observe that at this Reynolds number {Re ~ 6200) 
three-dimensional time-averaged azimuthal velocity is very close to the axisymmetric 
one both being upto 15% lower than the experimental data. The difference in Reynolds 
number is expected to play only a minor role in this discrepancy. 

The best fit of our (axisymmetric) computations (dash-dot line) with experimental 
data was realized when the boundary conditions (9) were perturbed with uniform 
random noise. The amplitude of the noise was 5% relative to the corresponding 
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maximums of axisymmetric solution without noise. The noise perturbation was applied 
for the part of the computational domain boundary of one spectral element long such 
that R± < r < R\+L g where L g = 0.0200 which is less than twice of the non-dimensional 
width of the gap between the inner cylinder and inner ring in the Princeton experiment 
equal to 0.0114 or 1.5 mm (Schartman 2008). 

The rational behind the noise perturbation of the boundary conditions was an 
attempt to model an effect of the centrifugally unstable flow in the gap between inner 
cylinder and inner ring. By an accident, boundary conditions on the inner cylinder 
surface (r = i?x) were also perturbed in this computation which turned out to be the 
best fit with experimental data. The effect of this perturbation of the inner cylinder 
surface boundary condition may be similar to a random blowing. This leads us to 
believe that a combination of the effects due to a run-out of inner cylinder and due 
to the centrifugally unstable flow in the gaps between the inner cylinder and inner ring 
may explain the discrepancy between the simulation and experiment (see also Schartman 
2008). 

The further work on comparison of the simulation with experiment is ongoing, and 
addition effort is needed to sort out the effects of run-out, centrifugally unstable flow 
in the gaps between cylinders and rings, vertical misalignment of horizontal surfaces of 
the rings, etc. 

3. 3. Torque and Angular Momentum Transport 

Also we have studied carefully the torque behaviour and associated angular momentum 
transport in the hydrodynamical setup of Princeton MRI liquid gallium experiment 
(Schartman 2008) as a baseline cases for our study of magneto- rotational instability 
(MRI) and MRI-driven turbulence. In order for the MRI experimental results to have 
a clear interpretation, the negative effects of the EC has to be minimized, and the 
torque amplification over the CCF torque Tq (15) with the increase of magnetic field 
can be linked directly to MRI enhancement of angular momentum transport (AMT). 
Therefore, the understanding of the torque behaviour and AMT in the baseline cases of 
hydrodynamical flow can not be overstated. 

An application of torque to the inner cylinder results in a flow that transports the 
angular momentum outward and attempts to reach a shear-free solid body rotation with 
a constant angular velocity. If the angular velocities of other boundaries are different 
form that of the inner cylinder, the resulting shear has to be maintained by application 
of torques to the boundaries in order to keep the rotation rates steady. In the context 
of MRI study, the primary interest is in the transport of the angular momentum from 
the inner to outer cylinder in centrifugally stable regime and, as shown below, in the 
minimization of EC and hence, in the reduction of the net contribution of torques exerted 
on the horizontal boundaries. Since the sum of all torques applied to the boundaries 
reflects the time increase of interior angular momentum, this contribution from the 
horizontal boundaries is equal to the sum of torques applied to the cylinders in a steady 
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Figure 6: The magnitudes of normalized torques applied to inner cylinder (Tj = T\/Tc, 
open symbols) and outer cylinder (T Q = —T 2 /Tc, filled symbols) versus Reynolds number 
for the cases with rings ( — Q— or ♦) and with lids (• -A- • • or A) in axisymmetric cases 
while the three-dimensional data are plotted with large symbols. Being the sum of 
magnitudes of torques applied to the inner and outer cylinder, the net torque exerted on 
the horizontal boundaries in the case of rings is significantly less than that in the case 
of lids making the former closer approximation to the CCF for which the net torque on 
horizontal boundaries is zero. 



state and has to be minimized for the successful MRI experiment. Note that in a case 
of unsteady flow, the time averaged torques are used when the statistically steady state 
is reached. 

The figure 6 shows the Reynolds number dependence of magnitudes of steady /time- 
averaged torque relative to the ideal CCF torque T~c (15) that has to be applied on inner 
and outer cylinders, Tj = ^ (open symbols) and T Q = (filled symbols), respectively, 
in order to keep constant boundary angular velocities (table 1). Being equal to unity for 
the case with periodic boundary conditions (CCF), the torque magnitudes are shown for 
the cases with rings (open diamonds with solid lines and filled diamonds) and lids (open 
triangles with dotted lines and filled triangles). All results are obtained in axisymmetric 
computations except for the data plotted with large open symbols that shows the results 
of fully 3D computations. Note that being zero in the ideal CCF, the difference between 
the torques applied to the inner and outer cylinder shown by open and closed symbols, 
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(a) Lids (b) Rings 

Figure 7: Steady state contour lines of effective angular momentum flux function \l> 
for the case of Re = 620 with lids in the range from -3.30 to 3.30 in the increment of 
0.31 (a) and with rings in the range from -1.35 to 1.35 in the increment of 0.21 (b). In 
the case of lids, most of the flux lines that originate from inner cylinder terminate at 
horizontal boundaries as opposed to the case of rings where they end up mostly at the 
outer cylinder which is similar to the CCF angular momentum transport between the 
cylinders. 

respectively, corresponds to the sum of torques exerted on the fluid next to the horizontal 
boundaries 

T t -T = (T 1 + T 2 )/T c 

due to zero net torque in steady/statistically steady state. Evidently, the setup 
with rings has an advantage of a smaller contribution to the net torque from the 
horizontal boundaries and of a smaller difference between inner and outer cylinder torque 
magnitudes over the setup with lids where EC is undisturbed. We also observe that in 
the range of Reynolds numbers considered the flow makes a transition from steady 
axisymmetric solution at Re = 620 to the unsteady one at Re = 6200 with small three- 
dimensional effects. Being more significant in the case with lids, three-dimensionality is 
expected to play an increasing role with the further increase of Reynolds number. 

In order to illustrate a spacial variations of AMT, we have computed an 
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effective angular momentum flux function defined in Appendix B by analogy with a 
streamfunction. The contours of the effective flux function shows the (flux) lines along 
which the angular momentum is transported, and the difference between the values of 
the flux function at two points gives the total flux across the segment of conical or 
cylindrical surfaces on which these points lie. Figure 7 shows steady state contour lines 
of constant increment for effective angular momentum flux function \l/ (B.12) for the 
case of Re = 620 with lids (a) and rings (b). For comparison, we note that the flux lines 
of (purely viscous) AMT for the ideal CCF (12) are the straight lines from inner to outer 
cylinder along z =const. Despite the fact that in both cases only a single line of 4> = 
(i.e. line of symmetry) is the same as in CCF case, the case with rings exhibits the similar 
transport of angular momentum along the flux lines that mostly originate at the inner 
cylinder and terminate at the outer cylinder in contract to the termination of the flux 
lines at the lids. The latter indicates that in the cases with lids, the angular momentum 
transport is mostly between the inner cylinder and the horizontal boundaries contrary 
to more desirable CCF-like transport between the cylinders observed in the cases with 
rings. Also note that the similarity between the shape of the flux lines away from the 
boundaries in figure 7 and the shape of vorticity contour lines and poloidal vector lines 
in figure 4 can be explained through creation of strong V r and V z components of velocity 
due to Ekman flows which affects AMT flux through advective contributions F® z and 
F% z , respectively, given by relations (B.5) and (B.6). 

In summary, if the ultimate objective is to achieve the flow with AMT as close to 
the ideal CCF as possible, the design with rings seems to have an advantage over the 
setup with lids. 

4. Conclusion and Future Work 

In this paper we have presented axisymmetric and fully three-dimensional Navier- 
Stokes calculations of circular Couette flow (CCF) in a cylindrical annulus as the first 
step in our study of magneto-rotational instability (MRI) and MRI-driven turbulence. 
Inspired by Princeton MRI liquid gallium experiment, we have computed the flow 
field in their experimental setup for realistic horizontal boundary conditions of 'lids' 
and 'rings' with the increase of Reynolds number through the onset of unsteadiness 
and three-dimensionality. The presented analysis of the flow field and angular 
momentum transport (AMT) allowed us to propose an explanation of the mechanism 
that determines the fate of the boundary flows and Ekman circulation (EC) as a result 
of a competition between the effects of 'centrifugal' rotation and pressure gradient set 
by rotation of, respectively, horizontal surfaces and bulk of the flow. In particular, with 
the appropriate choice of rotation rates of the horizontal rings that control an angle 
at which the vertical jets are launched near the ring boundaries, EC can be greatly 
diminished and CCF-like flow can be restored being more appropriate for the further 
experimental studies of MRI saturation and enhanced AMT. In addition, our numerical 
results compare favourably with the experimental data with the maximum deviation 
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below 15% being considerably smaller in the cases with 'noisy' boundary conditions. 
The future work, therefore, should involve higher Reynolds number computations with 
even more detailed modelling of experimental geometry that includes among others the 
effects of run-out of the inner cylinder, finite gaps between the cylinders and rings, and 
vertical misalignment of horizontal surfaces of the rings. 
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Appendix A. Traditional Explanation of Ekman Circulation 

Here we would like to comment that the traditional explanation of Ekman flows and 
circulation involves a balance between Coriolis and viscous forces in Ekman layers along 
rotating stressed boundary as viewed from a uniformly rotating reference frame (see 
?, ?). As we show below, this balance of Coriolis and viscous forces holds, for instance, 
in the case of a bulk flow outside the Ekman layers that is close to the state of solid 
body rotation at, say, angular velocity of fl 2 , when the centripetal pressure gradient and 
centrifugal forces cancel each other: 

p ~ -ny rnj + ^-o (A.i) 

2 2 2 dr 

(to the order of Ekman number E = A ^ L z ~ see ?) and Coriolis force is solely 
responsibly for the Ekman flow. In other flows like those considered in this paper, 
a contribution of the 'centrifugal' rotation to the Ekman flow does not reduce solely 
to Coriolis force but also includes the effects of centrifugal forces. Indeed, for the flow 
description in the non-inertial reference frame that rotates with constant angular velocity 
of the noslip boundary of the 'lids' f2 2 = ^2 e z , the centrifugal and Coriolis body forces 
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has to be added to the right-hand side of equations (3-5), 
/ = - fi 2 x (f2 2 x r) - 2 fl 2 x u 

f r = n 2 2 r + 2 Vl 2 u (A.2) 

along with the replacement of the inertial frame velocity V with the velocity in the non- 
inertial rotating reference frame u leaing to the following radial momentum equation in, 
e.g., case of axisymmetry (J^ = 0) 



du r _ u% ^ 2 1 r. u r ~\ dp 

_^l + u .V--^ = ft 2 y + 2ft 2 ue + — A (r:Z) u r -— -/ 

or 



st ^n 2!!2 " 6+ fJ-^ + rf + - (A - 3) 

Let us now compare this relation with equation (16) where we expand the 'centrifugal' 
term fl 2 r in powers of azimuthal velocity deviation from the state of solid body rotation 
of the lids, U e = V - rtt 2 : 

y r = Yl = (^ + U°) 2 =nlr + ( 20tUt + m (A.4) 



r r \ r 

Therefore, equation of radial momentum balance (16) rewritten with the expansion (A.4) 
coincides exactly with equation (A.3) in the case of axisymmetric flow when the 
difference between V r and u r and between Ue and u@ disappears due to the relationship 
between velocities in the inertial and non-inertial reference frames: 

V(r, 9, z,t) = r fi 2 e© + u(r, = 9 + Q 2 t, z, t) 

or 

V(r,z,t) =rtt 2 e e +u(r,z,t) (A.5) 

In a particular case of flows that are close to solid body rotation when considered in 
this description (A.3) equivalent to our earlier framework (16), the contribution of the 
centrifugal and pressure gradient forces cancels out (A.l) resulting in the steady state 
balance of Coriolis and viscous forces, 

1 d 2 u 

= 2n 2 u e + —^i + --- (A.6) 
Re ar z 

u 2 

where we omitted higher order terms including — 2 Vl 2 uq away from the axis of 
rotation r = due to vanishing u@ at the noslip boundary. In more general case, the 
additional contributions of the centrifugal forces Q 2 r and — to the Ekman flows has to 
be considered being a part of the 'centrifugal' term Q 2 r (16, A.3) that is balanced by 
centripetal pressure gradient and viscous forces in the Ekman flows along the lids that 
drive EC in the cylindrical annulus. 

Appendix B. Angular Momentum Flux and Flux Function 

Let us take the axisymmetric version of azimuthal momentum equation (4) with viscosity 
v instead of and rewrite it in terms of conservation of axial angular momentum 
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C = rVg. The summation of the axisymmetric versions of equations (4) and (6) 
multiplied by factors — r and rVg, correspondingly, gives 



-V, 



&Ve_ 
dr 



dVg We 

dz r 



+ v 



d 2 V e ldV e d 2 V 
H H 

dr 2 r dr dz 2 



Ve 



+ rV e 



dV r 8V, 



dr 



+ 



dz r 



dV ft dV ft 
rV r -^+rV z ^ + V r V e -v 



+ 



r dr 
dV r 



dz 

- rV e ^ + V r V e 



dr dz 
d[rV r V e ] , d[rV e V z 



+ 



d_ 

dz 

d_ 

dz 

— v 



dr dz 
dV e 



+ V r Ve - v 



d 2 V 6 
dr 2 



d_ 

dr 



v r 



dr 



\_dVe_ 
r dr 



dVe 
dr 



Vg 

rr.2 



v r- 



d^Ve 
dz 2 



v r 



d_ [V e 

dr 



-vr- 



dz 



tVqV z — vr 



dVe 



d_ 

dr 



dz 

(Va 



+ 



[rV r V e - 



vr 



2 JL (Ye)] b 



d_ 

dz 
ld_ 

r dr 



rVeV, — vr 



dr \ r 
dV e 



r'- - +V -V e 



dz ^ 
rV r Vg — vr 2 



| [rV r Ve-vr 2 j- r {^)] | d 
r dr 



rV r Vg — vr' 



d_ (V, 

dr 



l_ fV 

dr 



+ 



d_ 

dz 



rVgV z — vr 



dV 

dz 



(B.l) 



The above equation (B.l) reflects the conservation of axial or ^-component of 
angular momentum and can be condensed to 

dF zz 



or 



d Id 



|(r x V) + V,,„ • F 



dz 







= 



(B.2) 



(B.3) 



where total angular momentum flux tensor F and its components F rz and F zz are given 
for axisymmetric flows by 

(B.4) 
(B.5) 
(B.6) 

(B.7) 
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dVa 



vr 



dz 



-TT, 



(B.i 



where superscript a and v denote an advective and viscous contributions to the total 
flux of angular momentum, and r r9 and r zd are components of the stress tensor r. 



The Influence of Horizontal Boundaries on EC and AMT in a Cylindrical Annulus 23 



By analogy with poloidal streamfunction ip that satisfies axisymmetric version of 
continuity equation (6) due to, e.g., a definition 

* _i m 

an angular momentum flux function \l/(r, z) for steady flows can be introduced to satisfy 
the steady version of equation (B.2): 

As the lines of constant value of poloidal streamfunction ip are lines of constant poloidal 
flow rate with poloidal velocity vector being tangent to these streamlines, the lines of 
constant total angular momentum flux function \I/ are the flux lines along which angular 
momentum is transported, and the difference between these values at any two points 
gives the total flux of angular momentum transferred across the line that joints these 
points. Note that in the case of unsteady flow, the time-averaged quantities can be used 
in statistically steady state instead of instantaneous ones. Also note that similar to the 
streamfunction, the flux function W as a solution of equations (B.10) is defined up to a 
constant. We have used the point of half height on inner cylinder as a zero value for the 
flux function 

*(r = R 1 ,z = j y j = (B.ll) 

and for convenience, we define effective flux function $ as a flux function \1/ multiplied 
by circumference 2irr and scaled by the ideal CCF torque T c (15) 

27TT 

m(r,z) = —m(r )Z ) (B.12) 
T-c 

Finally, due to multiplication of flux function \1/ by r, the difference in effective flux 
function ^ between two points gives the total flux of angular momentum across the 
conical or cylindrical surface on which these points lie. Thus, in the case of e.g. inner 
and outer cylinder, it can be shown that torques T\ and T 2 scaled by the ideal CCF 
torque Tc are given by 

= H) - 0) ^ = V(R 2 , H) - 0) (B.13) 

Similar expressions hold for the torques applied to the horizontal surfaces. 



